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A new Quasiparticle Random Phase Approximation approach is presented. 

o _ 

. The corresponding ground state is variationally determined and exhibits a 

^ . minimum energy. New solutions for the ground state, some with sponta- 

. neously broken symmetry, of a solvable Hamiltonian are found. A non- 



iterative procedure to solve the non-linear QRPA equations is used and thus 

all possible solutions are found. These are compared with the exact results 

as well as with the solutions provided by other approaches. 

PACS number(s):21.60.Jz,21.60.Fw,21.60.Ev,23.40.Hc 
The random phase approximation (RPA) and its quasiparticle (qp) generalization 
(QRPA) are powerful tools for describing the collective degrees of freedom of many-fermion 
systems in various branches of physics like nuclear physics, solid state, plasma physics, 
atomic clusters. 

Usually the RPA and QRPA are presented by making use of boson expansion techniques. 
If the one body transition operator is expressed linearly in bosons, the many body Hamil- 
tonian with the two body interaction included becomes a quadratic polynomial in bosons 



describing harmonic motion of the selected degree of freedom. For the evident reasons the 
resulting approach is called a quasiboson approximation (QBA). The main drawback of this 
approximation scheme is that the QRPA (henceforth called as the standard QRPA) exhibits 
a collapse for a critical value of the strength of the attractive particle-particle force, i.e. the 
lowest eigenvalue becomes imaginary. In the case of the many body Hamiltonians used for 
the description of the transition rates of single and double beta decays the realistic proton- 
neutron force is expected to be close to its critical value. Due to this feature many theoretical 
works have been devoted to improving the standard QRPA by removing the ground state 
instabilities 

It is a common practice to investigate the validity of different approximation schemes 
within exactly solvable models, which describe the gross properties of the fermion many-body 
systems pl-p|,|7|,|8|, p!^ , p!5| . Recently, it has been shown that the QRPA does not develop any 



collapse and is in good agreement with exact results, if the Pauli exclusion principle (PEP) 
is consistently implemented in this approach . However, there is another group of studies 
saying that the collapse of the QRPA indicates a phase transition, i.e., a rearrangement of 
the nuclear ground state 0,^,^|T3|. In Ref. 0, in order to avoid the QRPA collapse, a new 
static ground state is defined by means of a semiclassical approach. A new collective mode, 
for a many-body system with proton-neutron pairing interaction, has been found beyond the 
point where the QRPA breaks down. This new solution requires the restoration of the PEP, 
objective that can be touched through a boson expansion technique. We note that the main 
difference between the semiclassical and the QRPA solutions consists in the description of 
the ground state. 

In this letter we point out new features for a many fermion system by improving the 
QRPA description of the ground state. A new ansatz for the QRPA wave functions is 
proposed, which allows us to minimize the QRPA ground state energy and, at the same 
time, to diagonalize the QRPA equations. The method proposed will be conventionally 
called as QRPA with an optimal ground state. One hopes that the present approach is a 
suitable tool for a realistic treatment of many body systems. It is worth mentioning that 



the correct description of the ground state of a many-body system is a steady interesting 
subject which has attracted the efforts of many groups P,^ p^p!3| . 

For the sake of simphcity we consider the proton-neutron monopole Lipkin Hamihonian 
which is exactly solvable 



Hp = e{Np + iV„) + XiA^A + X2{A^A^ + AA), (1) 

where Np (Nn) and A'^ are the proton (neutron) number and proton-neutron pair qp op- 
erators, respectively |^. The spaces of single particle states associated to the proton and 
neutron systems, respectively, are restricted to a single j-shell. The model Hamiltonian can 
be obtained from a Hamiltonian which in the particle representation, consists of pairing in- 
teraction for alike nucleons, and a monopole-monopole proton-neutron two body interaction 
of particle-hole and particle-particle types 0J|,|l^,|I|,|] . For the single j-shell, considered 
here, one obtains: 

where Up, Vp and m„, f„ are the coefficients of the Bogoliubov-Valatin transformation for pro- 
tons and neutrons, respectively. The parameters x and k are the strengths of particle-hole 
and particle-particle proton-neutron interactions, respectively. We consider the boson map- 
ping of the model Hamiltonian, following the Marumori recipe, i.e. the boson representation 
Hb is chosen so that its matrix elements (m.e.) between boson states are equal to the m.e. 
of Hp between the corresponding many fermion states The final result is [p|, [TO|] : 



Hb = auB+B + ao2(5+fi+ + BB) + a22B+B+BB 

+ ai3{B+BBB + B+B+B+B), (3) 



with 



an = {2e + Ai), ao2 = A2W (1 - ^), ^22 = 



a.,^A,[J(l-±)(l-i)-J(l-±)]. (4) 



The operator B"^ is a creation boson operator, while B is its hermitian conjugate operator. 

The QRPA treatment of Hp is determined by a RPA-hke procedure apphed to Hb- This 
consists in determining the phonon operator 

= X{B^ + t)-Y{B + t*), Q = {Q^y (5) 

as well as the vacuum state \rpa) and the one phonon state \Q) 

\Q>=Q^\rpa>, Q\rpa >= 0. (6) 

From the set of solutions depending on the parameters {t, t*) one depicts that one for which 
the state \rpa) has a minimum energy. In contradistinction to the usual QRPA, here the 
phonon operator involves a C number which results in having a nonvanishing expected value 
for Q, in the state \rpa), i.e. a "deformed ground state". 

The ansatz for the QRPA ground state \rpa >, defined as vacuum state for Q [see Eq. 
(|)], is taken as follows: 

|rpa>=e**^-*^'e^^'^'-^*^^|>. (7) 

Here | > denotes the vacuum state for the boson operator B, which might be the mapping 
of the uncorrelated BCS ground state, t, z are complex parameters: t = re*^, z = pe^'^. The 
ansatz (|^ obeys the second eq. (|^) provided the following equation holds: 

Y 

T[r, p, 0) = sinh (2p) - —e''^ cosh (2p) = 0, (8) 

A 

where the function depends implicitly on r and 6 variables, by means of X and Y ampli- 
tudes. It is worth to note that the functions describing the limiting cases (p, 0) = (0, 0) and 
(r, 6) = (0, 0) were earlier considered as trial states in a time dependent formalism: 

\rpa) = e-'^e-*^', \rpa >= ^ e'^^'^'l > (9) 

-y/cosh (2p) 

with d = e*"^ tanh (2p)/2. Due to this feature one expects that the present trial function 
accounts for new correlations in the ground state. The phonon amplitudes X, Y satisfy the 
QRPA equations 



AX + BY^ ErpaUX, 

B*X + AY = -ErpaUY, (10) 

where the RPA energy Erpa is the excitation energy for the RPA state \Q), i.e. E^pa — 
El — Eg s.- The RPA matrices A,B,U, are determined in the usual manner and have the 
expressions: 

A = < rpa\ [B, Hb, B^] \rpa > 

= Oil + 4022^^ + Gaisuv cos 4> + r^[4a22 + 6013 cos (26*)], (11) 
B = < rpa\[B,HB,B]\rpa > 

= 2ao2 + 2a22Uv[cos (j) + isin0] + Gai^v'^ + r^(2a22[cos (26*) + isin (26*)] + 6013), 
U=< rpa\ [B, B^] \rpa >= 1. (12) 

The double commutators are defined as 2[A,B,C] = [A, [S,C]] + [[A, S],C]. Also, the 
notations u — cosh (2p) and v — sinh (2p) have been newly introduced. The matrix A is 
real, while B might be a complex number. In this letter we consider only situations where 
the solutions for the X, Y amplitudes are real. The arguments which support our choice 
is that the resulting phonon state should be the image, through the boson mapping, of a 
state in the fermionic space. However the amplitudes of the eigenstates of Hp, in the basis 
{(^4''^)"}, are real numbers. Thus the possible solutions for 9 and are 9 — n7r/2 and (f) — nn 
(n — 1,2,...). The additional parameters r and 9 entering the QRPA eigenvalue problem will 
be fixed by requiring that the expectation value of Hb in the RPA ground state reaches its 
minimum value Eg s., in the space of these parameters. One finds 

H{r) = < rpa\HB\rpa >= Co + Csr^ + (13) 
Co — aiiv^ + 2aQ2UV cos + a22[u^v'^ + 2i'^] 

+6ai3uv^ cos (f), 
C2 = an + 2ao2 cos {26) + 022 [4f^ + 2uv cos (261 + (p)] 

+6ai3 \uv cos 4> + v'^ cos (2^)] , 

5 



C4 = a22 + 2ai3Cos(2^). (14) 
Several minima for Eg ^ are to be mentioned: 

a) r = and 9 = anything. This solution corresponds to the RPA ground state given by 
the second Eq. (^). 

b) r = y'-C2/2C4 , ^ = 0, c) r = ^^C^/2Cl, 9 = 7i/2, 
d)r = ^-C2/2Ca, 9 ^ n-K/2 (n=l,2,3 ...)• 

The case d) yields a complex B and therefore is not considered here. The solutions b), c) 
and d) are specific for the ansatz (|^). If only the quadratic terms in bosons are considered 
in Hb (i.e., 022,^13 = 0) then only the solution a) survives. The equations (10), (8) and 
one of a),b),c) are to be simultaneously solved. This non-linear problem is usually treated 
iteratively. However, this procedure does not guarantees that all solutions of this system of 
equations, are found. The procedure used in the present paper is as follows. First, we insert 
the solution for r (r = or = — C2/2C4) in the expressions defining the matrices A and B 
by Eqs.(ll), (|12]). Determine the ratio Y/ X from any of equations (|T0|) and insert the result 
in Eq. (|^). Then Eq. (^ provides p which is to be inserted in the equations for the RPA 
matrices. With the matrices A, B fully known one determines the QRPA eigenvalue and 
eigenvectors, by solving the equations (10). This procedure is repeated for each minimum 
of Eg s. [the cases a), b) and c) mentioned above]. We note that solutions with (p, (p) and 
(— p, ± tt) are degenerate. 

The numerical application is performed for a system of 4 protons and 6 neutrons moving 
in a j = I shell. The strength parameters defining the model Hamiltonian are e = IMeV 
while the two body interaction strengths k and x ^i^re re-scaled as in Refs. k, ^ k' = 

2Q K, X ^ x' = 2^^ X- We adopt x' = 0.5 while k' is considered as a free parameter in 
the interval (0, 3)MeV. In Table | we present the results for the QRPA and ground state 
energies corresponding to three values of k'. If n' = 0.5 there are only two solutions of 
type a), associated with the ansatz (^ (second equation). The solution a2) is similar to 
that obtained in |I| within the EPP QRPA formalism. The solution al), with small QRPA 
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energy and a very large value of Eg ^., has not been seen before. Its finding is a merit of 
the adopted numerical procedure. For k' = 1.5, 2.5 one finds solutions of type c) which are 
specific for the new ansatz of the RPA wave functions. 

In Fig. we present the potential energy surface as function of r for the solutions a2), 
cl), c2) and c3), assuming k' = 2.5. Note that for the solution c) H{r) has two minima, at 
r = ±^-C2/(2C4). 

The meaning of the negative r consists of that the result for energy is not affected if we 
change the ansatz by replacing t with —t. 

In Fig. H the values of p corresponding to different solutions are presented. Within 
the QBA, i.e. the quadratic bosonic Hamiltonian is considered, the value of p is increasing 
rapidly in the vicinity of the QRPA collapse. The solution a2) does exist within the whole 
interval of k'. This is caused by that the PEP principle included in the higher order bosonic 
expansion Hamiltonian prevents the collapse of the QRPA The solution al) was found 
within the range 0.0 < n' < 1.29 and continues further as cl) solution. The solutions c2) 
and c3) start with k' = 1.53 where they have equal values. For large value of k' the solution 
c2) is close to a2). It is worthwhile to notice that in the case of c3), with the lowest value 
of Eg s and large k' (see Table |), the corresponding values of p are rather small. One may 
say that the corellations in the ground state are mainly induced by the "deformation" of 
the system, i.e. the stationary value of r. For n' > 2, p is slightly negative. This branch is 
obtained by shifting the phase with vr which is equivalent to changing the sign of p but 
keeping unchanged. 

In Fig. 1^ we plotted the ground state energy Eg g, the QRPA energy Erpa and the first 
excitation energy Ei obtained by the standard QRPA, QRPA with an optimal ground state 
and by diagonalizing Hp, as function of k'. Only solutions with low ground state energy 
are drawn. One remarks that by including the higher order boson terms in the Hamiltonian 
Hb and considering the standard ansatz for \rpa) (see the second Eq. (9)) [solution a2)], 
there is no collapse of the QRPA solution. This solution coincides well with that provided 
by the EPP QRPA formalism |]l|. However, for large values of n' {n' > 1.6) this solution 
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approximates poorly the exact result for Eg,s ■ On other hand the new solution c3) is in 
an excellent agreement with the exact Eg g for this range of k'. In fact, the exact result 
for Eg,s is very well described by matching the approaches a2) and c3). Despite the fact 
the c3) approach simulates excellently the exact result for the ground state energy, the two 
treatments predict QRPA energies which deviate from each other by a large amount. 

In summary, a new ansatz for the QRPA wave functions was proposed. The additional 
new factor in this ansatz depends on a complex parameter t which is fixed by requiring 
that the corresponding expectation value of the model Hamiltonian Hb is minimum. The 
present paper constitute the first attempt in the literature, to determine the QRPA ground 
state (the vacuum state for the Q operator) as a solution of a variational equation. The 
numerical application shows clearly new effects caused by the presence of higher order boson 
terms in Hb and by the complex structure of the ansatz for the state \rpa). Both features 
mentioned above prevent the QRPA to collapse. Several solutions are found in the interval 
of the interaction strength k' beyond the critical value. Each solution corresponds to a 
certain type of minimum in the (r, 6) variable and by that to a certain symmetry of the 
wave function. One may say that the present approach accounts, in an unified fashion, for 
the complementary features described by the EPP QRPA and the semiclassical formalism, 
which defines a new ground state for those k' where the standard QRPA fails. 
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TABLES 

TABLE I. The predicted QRPA (Erpa) and ground state {Eg,s) energies for k' = 0.5, 1.5 and 
2.5, x' — 0-5 and j = 9/2. The parameters determining the RPA wave functions and the ground 
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FIG. 1. Expectation value of the Hamiltonian Hb as a function of the parameter r. The 

solutions labeled by a) and c) are those presented in the text below Eq. (14). The adopted values 
of x' and k' are 0.5 and 2.5, respectively. 
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FIG. 3. Energies of the ground state Eg,s, of the first excited state Ei and the excitation 
energies Erpa{= E\ — Eg,s) provided by diagonaUzing Hp, by the standard QRPA and by the 
QRPA with an optimal ground state [ al), c2) and c3)], respectively, are plotted as function of k'. 
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